Unconventional pairing originating from disconnected Fermi surfaces in 

superconducting LaFeAsOi-^Fj, 
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For a newly discovered iron-based high T c superconductor LaFeAsOi-ajFa;, we have constructed 
a minimal model, where inclusion of all the five Fe d bands is found to be necessary. Random-phase 
approximation is applied to the model to investigate the origin of superconductivity. We conclude 
that the multiple spin fluctuation modes arising from the nesting across the disconnected Fermi 
surfaces realize an extended s-wave pairing, while d-wave pairing can also be another candidate. 
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Understanding the mechanism of unconventional su- 
perconductivity (SC) has been one of the most challeng- 
ing problems in condensed-matter physics. There is a 
renewed fascination with a recent discovery of SC in 
an iron-based superconductor LaFeAsOi-zFa;,!!! which 
is likely to provide a fresh avenue for such a challenge. 
LaFeAsO belongs to the family of quaternary oxypnic- 
tides LnMPnO (Ln=La, Pr; M=Mn, Fe, Co, and Ni; 
Pn=P, As), which was originally fabricated by Zimmer 
et al. and Quebe et a/.0, 0] For this family of com- 
pounds, Kamihara et al. first reported that LaFePO ex- 
hibits SC with T c ~ 3K, which was raised to T c ~ 7K by 
F doping. 0] SC has also been found in nickel-based com- 
pounds with the same structure. [5j Very recently, Kami- 
hara et al. have come up with the discovery of SC in 
LaFeAsOi-ajF^, where the F doping with x ~ 0.11 leads 
to a remarkable T c ~ 26K. 

The high value of T c itself, confirmed also by Chen et 
aL,0| suggests a possibility of unconventional SC, but 
direct evidences are accumulating: A specific heat mea- 
surement in magnetic fields shows that the coefficient 
7 displays a VH behavior. p| A point-contact conduc- 
tance measurement exhibits spectra with a distinct zero- 
bias peak,Q suggestive of the presence of sign change 
in the gap function, [jj [n| [TJ, [lj| The starting ma- 
terial, LaFeAsO, is a bad metal with some anomaly in 
the resistivity around 100K. 1] As the system becomes 
metallic upon F doping, the uniform susceptibility ex- 
hibits a Curie- Weiss behavior. Anomalies in the normal- 
state transport properties have also been reported for 
doped systems. [1 3} 

Theoretically, first-principles band structure has been 
obtained for LaFePO 14 L and more recently for LaFeAsO 



correlation enhances the crystal field splitting, which 
leads to a band-semiconducting behavior in accord with 
the experiment. [l6j Local spin-density calculations for 
LaFeAsO show that the system is around the border be- 
tween magnetic and nonmagnetic states, with a tendency 



and related materials ljj H, 12, EH- These band struc 
tures are metallic with five pieces (sheets) of the Fermi 
surface in the undoped system, which contradicts with 
the experiment for the undoped LaFeAsO. [l| However, 
a dynamical mean-field study shows that the electron 



toward fcrromagnctism and antiferromagnetism. la, LL7| 
It is also pointed out that the electron-phonon coupling 
is too weak to account for T c = 26K. 0, E3| 

Given this background, the purpose of the present Let- 
ter is to first construct a microscopic electronic model for 
LaFeAsOi-zF-r, which then serves as the basis for iden- 
tifying the possible mechanisms why this material favors 
high-T c . The minimal model has turned out to contain 
all the five Fe d orbitals, to which we have applied the 
random-phase approximation (RPA) to solve the Eliash- 
berg equation. We shall conclude that a peculiar Fermi 
surface consisting of multiple pockets and ensuing mul- 
tiple spin-fluctuation modes realize an unconventional s- 
wave pairing, while d-wave pairing can also be another 
candidate. 

LaFeAsO has a tetragonal layered structure, in which 
Fe atoms are arrayed on a square lattice. Due to the 
tetrahedral coordination of As, there are two Fe atoms 
per unit cell. Each Fe layer is then sandwiched between 
LaO layers. The experimentally determined lattice con- 
stants are a = 4.03552Aand c = 8.7393A, with two inter- 
nal coordinates zj, a — 0.1415 and Za s — 0.6512. We have 
obtained the band structure (FigOJa) inset) with the 
Quantum-ESPRESSO package [ljj, and then construct 
the maximally localized Wannier functions (MLWFs) [2(| • 
These MLWFs, centered at the two Fe sites in the unit 
cell, have five orbital symmetries (orbital l:d 3 z 2 -R 2 , 
2-dxZi 3'dyz, ^'-dx 2 -Y 2 , 5'-dxY, where X,Y,Z refer to 
those for the original unit cell). We can note that the 
two Wannier orbitals in each unit cell are equivalent in 
that each Fe atom has the same local arrangement of 
other atoms. We can thus take a unit cell that con- 
tains only one orbital per symmetry by unfolding the 
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TABLE I: Hopping integrals t( Ax, Ay; fj,, v) in units of 
O.leV. [Aa;, Ay] denotes the in-plain hopping vector, and 
(fi, v) the orbitals. a y , I, and ad corresponds to 
t(Ax, -Ay;/i, u), t(-Ax,-Ay,fj,,v), t(Ay, Ax; fi, i>), respec- 
tively, where ' + ' and ' + (//, v')' in the row of (/j, v) mean that 
the corresponding hopping is equal to ±t(Ax,Ay; jj,,u) and 
±t(Ax, Ay; //, v 1 ), respectively. This table, combined with 
the relation t(Ax, Ay; fi, v) = t(~Ax, —Ay; gives all the 
in-plain hoppings > O.OleV up to fifth neighbors. 



Brillouin zone (BZ),[2J| and we end up with an effec- 
tive five-band model on a square lattice, where x and y 
axes are rotated by 45 degrees from X-Y (FigOJb) in- 
set), to which we refer for all the wave vectors hereafter. 
The in-plain hopping integrals t(Ax, Ay, Az = 0; ji, v) 
are displayed in table U where [Ax, Ay] is the hopping 
vector, and y,, v label the five Wannier orbitals. The on- 
site energies for the five orbitals are (ei, £2, £3, £4, £5) = 
(10.75, 10.96, 10.96, 11.12, 10.62) eV. With these effective 
hoppings and on-site energies, the in-plane tight-binding 
Hamiltonian is given in the form 

ij flu a 

+ t(xj - x l ,y j - yi, v, y)c] va c illiy + ^2 £ ^ n n^i (!) 
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where cj creates an electron with spin a on the /x-th 

orbital at site i, and — c\ c i/xa .. We define the 

band filling n as the number of electrons/number of sites 
(e.g., 71 = 10 for full filling). The doping level x in 
LaFeAsOi-zF^ is related to the band filling as n = 6 + x. 

In the obtained band structure in Fig[TJa), we no- 
tice that the five bands are heavily entangled, reflecting 
strong hybridization (see tabic [TJ) of the five 3d orbitals, 
which is physically due to the tetrahedral coordination 
of As atoms around Fe. Hence we conclude that the 
minimal electronic model requires all the five bands. In 
FiglUb), the Fermi surface for n = 6.1 (corresponding to 



x = 0.1) obtained by ignoring the inter-layer hoppings 
is shown in the two-dimensional unfolded BZ. The Fermi 
surface consists of four pieces (pockets in 2D): two con- 
centric hole pockets (denoted as a±, a-i) centered around 
(k x ,k y ) — (0,0), two electron pockets around (tt,0) 
or (0,7r) (/3a), respectively. a« (/%) corresponds to the 
Fermi surface around the TZ (MA) line (in the ori gina l 
BZ) in the first-principles band calculation. 3, IH 13] 
Having constructed the model, we now move on to the 
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FIG. 1: (color online) (a) The band structure of the five-band 
model in the unfolded BZ, where the inter-layer hoppings are 
included. To compare with the ten-band model (red lines in 
the inset; the symbols are the present LDA results), note the 
original (dashed lines) and the unfolded (solid) BZ shown in 
Fig.(b). (b) Fermi surface for n — 6.1 (with the inter-layer 
hoppings ignored) , with the arrows indicating the nesting vec- 
tors. Inset depicts the original (dashed) and reduced (solid) 
unit cell in real space. 
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RPA calculation. We again adopt the 2D model in which 
the inter-layer hoppings are neglected. [22j For the many- 
body part of the Hamiltonian, we consider the stan- 
dard interaction terms that comprise the intra-orbital 
Coulomb U, the inter-orbital Coulomb U', the Hund's 
coupling J, and the pair- hopping J'. All the calculation 
is done in the orbital representation. Details of the multi- 
orbital RPA calculation can be found in e.g. ref.[23, 24 1. 
The modification of the band structure due to the self- 
energy correction is not taken into account, on which we 
comment later. In the present case, the Green's func- 
tion is a 5 x 5 matrix, while the spin and the orbital 
susceptibilities become 25 x 25 matrices. The Green's 
function and the effective pairing interactions, obtained 
from the susceptibilities, are plugged into the linearized 
Eliashberg equation, and the gap function in a 5 x 5 ma- 
trix form along with the eigenvalue A is obtained. T c 
corresponds to the temperature where A reaches unity. 
32 x 32 fc-point meshes and 1024 Matsubara frequencies 
are taken. We find that the spin fluctuations dominate 
over orbital fluctuations as far as U > U' , so we focus on 
the spin susceptibility. We denote the largest eigenvalue 
of the spin susceptibility matrix for ioj n — as Xs(k). 
The gap function matrix at the lowest Matsubara fre- 
quency is transformed into the band representation by 
a unitary transformation, and its diagonal element for 
band i is denoted as <j>i(k). 

Let us first look at the result for \s for U = 1.2, 
U' = 0.9, J = J' = 0.15, and T = 0.02 (all in units 
of eV) in Fig[2^a). The spin susceptibility has peaks 
around (k x ,k y ) = (it, 0), (0, it) and also a ridge-like 
structure from (ir,ir/2) to (tt/2,tt). This in fact reflects 
the Fermi surface in FigQJb), where we have two kinds 
of nesting vector: ~ {%, 0), (0, 7r) across a and /3, and 
~ (-7T, 7r/2), (tt/2, 7r)across j3\ and /?2- A good nesting en- 
hances tendency towards magnetism. In the RPA (where 
the self-energy correction in the Green's function is ne- 
glected), we have to take U as small as 1.2eV to ensure 
magnetic ordering does not take place in the temperature 
range considered. 

For SC, we show in Fig[2]Jc)(d) the gap function for 
bands 3 and 4 (as counted from below), together with 
the Fermi surface of each band. At this temperature 
(T = 0.02) ,the eigenvalue of the Eliashberg equation is 
A = 0.96. [25j The gap is basically s-wave, but changes 
sign between the Fermi surface of band 3 (0.2) (and also 
band 2; oti, not shown) and those of band 4 (/3i , /?2) , 
namely, across the nesting vector ~ (tt, 0), (0, ir) (M point 
in the original BZ) at which the spin fluctuations de- 
velop. Such a sign change of the gap between inner 
hole and outer electron Fermi pockets is analogous to 
those in models studied by Bulut et aZ.,[26| and also 
by two of the present authors. 27. |28| It is also remi- 
niscent of the unconventional s-wave pairing mechanism 
for Na a; Co02 • 2/H2O 29] proposed by four of the present 
authors. (30| After completion of the present study, we 



have come to notice that a recent preprint by Mazin et al. 
also concludes an s-wave pairing in which the gap changes 
sign between a and (3 Fermi surfaces .[3lj| as schematically 
shown in the upper panel of Figl^b). For the present set 
of parameter values, in addition to this sign change, we 
find that the nodes of the gap intersect the (3 Fermi sur- 
face. This is because the spin fluctuations due to the 
01-13% nesting favor a sign change in the gap between (3\ 
and 02 Fermi surfaces. In fact, we have found that this 
nodal line moves out of the (3 Fermi surface for the pa- 
rameter values for which the spin fluctuations due to the 
(3\-(32 nesting become less effective, e.g., for U — U' . In 
that case, the gap becomes closer to the upper panel of 
FigHb).[33 

We have so far focused on the diagonal elements of the 
gap matrix in the band representation, but to be more 
precise, we have to consider the off-diagonal (interband) 
elements in order to make accurate comparison with ex- 
periments. The off-diagonal elements in the present case 
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FIG. 2: (color online) RPA result for the spin susceptibility \ s 
(a), the gap functions 4>3 (c) and <j>4 (d), y (0</>t) 44 (e) for U = 
1.2, U' = 0.9, J = J' = 0.15, n = 6.1 and T = 0.02 (in eV). In 
(c) and (d), the black (green) solid lines represent the Fermi 
surfaces (gap nodes). In (b), the fully gapped extended s- 
(upper panel) and d I 2_ !/ 2-wave gaps are schematically shown. 
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turn out to be not negligibly small due to the heavy en- 
tanglement of the bands. One way to look at this effect 

is to calculate the quantity J (<^t) 44) where <fi is the gap 
matrix and 44 denotes the diagonal element of band 4. 
As shown in Fig^e), we find that this quantity is finite 
over the entire BZ, but a remnant of the nodal lines of 
the diagonal element still appears as a valley that inter- 
sects the j3 Fermi surface. In this sense, we can say that 
the magnitude of the gap varies along the f3 Fermi sur- 
face (becomes large at points far from the BZ edge) if the 
spin fluctuations arising from a- (3 and /?i-/?2 interactions 
have competing strength. The degree of the variation 
of the gap in the actual materials may be determined 
experimentally from the density of states, e.g., tunneling 
spectroscopy, or directly by angle resolved photoemission 
studies. 

In the above, we mainly considered the possibility of 
unconventional s-wave pairing. On the other hand, if the 
a Fermi surfaces are absent (or less effective), the sim- 
plest form of the gap would be the d x 2 _ y 2 -wave pairing 
{dxY in the original BZ), where the gap changes sign 
between (3\ and (32 Fermi surfaces as shown in the lower 
panel of FigJ^b). To check this, we have performed an 
RPA calculation on (i) the present model with n = 6.3 
and (ii) a model where we artificially shift the crystal 
field splitting to let the a Fermi surfaces disappear for 
n = 6.1. In both cases, we indeed obtain the d x 2_ y 2- 
wave. Since the band structure generally changes from 
the LDA result due to correlation effects[16| or a band fill- 
ing different from the formally expected value, we leave, 
at the present stage, this d-wave state as another candi- 
date for the pairing symmetry in this material. 

Many other interesting problems remain for future 
studies. Spin fluctuations and SC should be studied by 
taking into account the self-energy correction, for which 
a fluctuation exchange study is underway[32] It is also 
intriguing to investigate whether the present unconven- 
tional gap can quantitatively account for the specific- 
heatQ and point-contact conductance [|| results. Also, 
further insight into the origin of the high T c SC in 
LaFeAsOi-aTa, may be obtained by performing similar 
microscopic studies on LaFePOi-aT^T] or LaNiPO.jH 
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